Metallic-to-insulating transition in disordered graphene monolayers 
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We show that when graphene monolayers are disordered, the conductance exhibits a metallic- 
to-insulating transition, which opens the door to new electronic devices. The transition can be 
observed by driving the density or Fermi energy through the mobility edge. At the Dirac point 
the system is localized, whereas at higher densities there is a region of metallic behavior before the 
system becomes insulating again at higher densities. The region of metallic behavior depends on 
the disorder strength and eventually vanishes at high disorder. This result is quite unexpected since 
in square lattices, scaling theory predicts that this metallic region does not exist in two dimensions, 
in contrast to graphene, where the lattice is a honeycomb. 



In most active electronic devices, the conductivity can 
be tuned from conducting to insulating by using a gate. 
With the recent discovery of graphene monolayers [![, 
and their potential for electronic devices Q, it is impor- 
tant for the conductivity to change substantially with the 
gate voltage. However, while in clean graphene monolay- 
ers, the conductivity changes almost linearly with gate 
voltage, the off or minimum conductivity is still rela- 
tively large (of the order of the unit of conductance e 2 /h) 

ass. 

Here we consider the situation of graphene nanorib- 
bons (with a honeycomb structure) of width W and 
length L contacted by two large normal metals, assumed 
to have a square lattice. This leads to a typical two- 
terminal configuration as shown in figure [TJ where the 
contacts are assumed to be perfect. The disorder poten- 
tial is taken to be only onsite and uncorrelated. This is 
identical to the large body of work on Anderson localiza- 
tion in tight binding models Q, which show a metal- 
insulator transition in dimensions strictly higher than 
2 0, 0, [l^, 11 1 . In two dimensional systems with a 
square lattice no metal-insulator transition exists for non- 
interacting electrons JH, H3 |. except if correlations in the 
disorder are present [131 ]. 

In graphene, with its honeycomb lattice, the band 
structure was first studied theoretically by Wallace [14 1 
using a tight binding Hamiltonian, where t ~ 3eV is the 
hopping energy between two nearest neighbors and V n 
the onsite energy. At the band edges, the band structure 
is very similar to the band edge of square lattices. At the 
center of band, however, and for V n = Wallace showed 
that there is a linear dispersion at the band center (Dirac 
point) at two points in the reciprocal space leading to a 
two-valley degeneracy. This has important implications 
on scattering properties, such as a suppression of intra- 
valley backscattering. The linear dispersion close to the 
Dirac point leads to dramatic new physics, such as an 
anomalous quantum hall effect and Dirac fermions [|| . It 
is the absence of a gap at the Dirac point, which causes 
the conductivity to not vanish at that point 0, H, 0], 
which is a potential roadblock for applications in active 
electronics. 
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FIG. 1: Top: Image of a graphene nanoribbon with two 
attached contacts. Main graph: Length dependence of 
— (log(G)), — 21og(G) and log(_R)/2 for different energies and 
widths. The average is taken over 200 configurations and 
V = 1.5. The solid lines represent linear fits of the data for 
L 2> Lc, the slope of which gives the inverse localization 
length Lc- The horizontal line illustrates the value of Lh- L 
and W are given in units of graphene atoms. 



Disorder plays an important role in graphene devices 
and can be due to ripples HI, 1|| , defects in the substrate 
17 1 and surface effects, such as partial hydrogenization 
181 ] . To reduce disorder scattering, suspended graphene 
devices [l|| have been considered as well as the use of 
other substrates [20| . 

Localization in graphene nanoribbons: Early numerical 
studies in disordered honeycomb lattices, were limited to 
the Dirac point and showed that the states are local- 
ized at this point [l2j]. These results were obtained by 
evaluating the localization length (Lc) from the small- 
est Lyapounov exponent of a finite width (W) ribbon 
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using a transfer matrix approach. Lc is then studied 



as a function of W. Using scaling arguments lo- 
calized states are identified when L c /W decreases with 
W, whereas extended states are characterized by Lc/W 
increasing with W . The point where there is no de- 
pendence is then inferred to as the critical point, where 
the localization-delocalization transition occurs. In two 
dimensional square lattices, only localized states were 
found for all energies 21 1. 

Here we use a very similar approach, but instead 
of considering a transfer matrix, we evaluate the two- 
terminal conductance of the system, which allows us 
to directly evaluate the transport properties of the sys- 
tem. For a given width of the graphene nanoribbon, 
we evaluate the zero-temperature two-terminal conduc- 
tance (G) of the disordered graphene attached to metallic 
(square lattice and non-disordered) leads using an itera- 
tive Green's function technique. Since the conductance 
depends on the given disorder configuration we consider 
a configurational average (•) by averaging over many dis- 
order configurations. We used two disorder distributions, 
either uniform with — V/2 < V n < V/2, or binary (V n = 
±V/a/12) both characterized by (V n V m ) = S n , m V 2 /12. 
The two distributions give identical results, in contrast 
to one dimension, where large distribution dependent en- 
semble averaged conductance fluctuations exist [22|. We 
assume the leads to be disorder free and infinitely long 
and much wider than the disordered graphene rectangle. 

In the presence of disorder, the conductance will even- 
tually vanish exponentially with the length of the sys- 
tem. This is illustrated in figure [TJ where we show 
— (log(G)) for different values of the Fermi energy in 
units of t. Because of the symmetric band structure 
around the Dirac point, all results are symmetric around 
E = 0. Two regimes can be identified: (i) the ballis- 
tic regime, when L <C Lc and (ii) the localized regime, 
where -(log(G)) ~ L/Lc for L > L C - In the bal- 
listic regime {L <C Lc), the conductance is dominated 
by mesoscopic conductance fluctuations [U 
SG < G, hence 



24|, where 



(log(G)) ~-log(G) ~log(i?). 



(1) 



R = l/G is the two terminal resistance of the device and 
SG is the standard deviation of G. In the localized regime 
(L ^> Lc), on the other hand, where the conductance 
vanishes exponentially with the length of the system, we 
have SG ^> G. Because of the statistical properties of 
the conductance of a quasi-one dimensional system, this 
yields 



(log(G)) ~ -21og(G) ~ log(i?)/2 ~ L/Lc - 



(2) 



a is a parameter, which is close to unity for G in units 
of 2e 2 /h as shown in figure [TJ Relation (2) becomes ex- 
act in the limit where L — > oo. Hence, the localization 
length can be extracted using any of the average trans- 
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FIG. 2: Left: The dependence on width of the ratio Lc/W 
for different values of the energy. Right: The dependence of 
the ratio Lc/W on energy for different values of the width. 
The energy Ec labels the critical energy at which the width 
dependence of the ratio changes sign, indicative of a metallic- 
to-insulating transition. 



port quantities ((log(G)), (G) or (R)), but the conver- 
gence of (log(G)) is much faster as shown in figure [JJ 

While equation (2) requires L ;» L c in order to ex- 
tract Lc accurately, an approximate Lc can also be ob- 
tained by looking at the crossover from ballistic to lo- 
calized, which will happen when equations (l) and (2) 
both hold, i.e., when each term is zero. This corresponds 
to SG ~ G, since SG is close to one in units of 2e 2 /h 
[23I ] . Hence, defining Lh as the length which minimizes 
min{(log(G(L))) 2 + (log(G(L))) 2 + (log(i?(L))) 2 } => L H , 
we obtain a length, which characterizes the crossover 
from ballistic to localized. It turns out that for all values 
of interest, L# ~ Lc within 10%. 

Scaling behavior: In order to determine if a state is 
exponentially localized in the two dimensional limit, we 
need to evaluate Lc/W as a function of W. Lc is always 
finite in the quasi one-dimensional limit, since all states 
are localized. For the two-dimensional case, Lc/W is the 
relevant quantity, since ultimately we are interested in 
the average conductivity a = (G)L/W for infinite width 
and length. By leaving the aspect ratio constant and 
equal to one, i.e., L — W when taking the limit to large 
sizes, we obtain a ~ e -w/2L c ^ smce (Qj ^ e -L/2L c f rom 

equation (2). Hence, if with increasing width Lc/W — > 
this leads to a — > and we have an insulator. On the 
other hand, if Lc/W — > 00, for increasing width, this 
implies that there are no exponentially localized states 
and we refer to this state as metallic. 

We now apply this analysis to our graphene device and 
show the results for Lc/W in figure O We observe that 
at the Dirac point, the ratio Lc/W monotonously de- 
creases with the width, which implies that the system 
is insulating at the Dirac point. This is consistent with 
earlier results [l2, 26|, 27 1. More interestingly, away from 
the Dirac point, the dependence of Lc/W is more com- 
plicated and the dependence becomes non-monotonous. 
Indeed, at small widths the ratio decreases at first, before 
increasing again at larger values of the width, which is a 
signature for metallic behavior. This increase of Lc/W 
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FIG. 3: Contour plots of the ratio L c /W as a function of the 
energy and width. Top graph is for V = 1.5, bottom left for 
V = 2 and bottom right for V = 2.2. 



only occurs in a small window of energies, between E = 
and E = Ec as shown in figure O This is typical for 
metal-insulator transitions as seen in three dimensions 
11]. However, this is quite unexpected in two dimensions 
without a magnetic field. Indeed, numerical studies on 
disordered systems with square lattices show no metallic 
behavior [12j, i.e., Lc/W always decreases with size. 

Interestingly, the continuous increase of Lc/W over 
a wide range of widths close to the Dirac point, which 
is a signature for a metallic behavior, is correlated with 
the existence of additional channels. Indeed, for a non- 
disordered system close to the Dirac point, there is only 
one quantized transport channel, leading to a quantized 
conductance of 2e 2 /h [28| if the width of the nanorib- 
bon is smaller than W\u = 2ttHvf/ Ep, where vp is the 
Fermi velocity and Ep the Fermi energy. The value of 
Wm is shown in figure [2] for Ep — — O.lt and is corre- 
lated to a jump in L c /W . This is quite different from the 
square lattice case, where the addition of another chan- 
nel is always correlated with a trough of Lc [S] and 
suggests the two situation to be very different in nature 
and indicates that the reason for this delocalization is 
intimately related to the linear dispersion, which leads 
to a suppression of intra-valley back-scattering. In terms 
of the beta function, Nomura and coworkers argued for 
the topological delocalization of two-dimensional mass- 
less Dirac fermions |30j |. 

A more detailed analysis of the dependence of Lc/W 
is provided by the contour plots shown in figure [3l For 
V = 1.5 the phase space in energy and width where 
Lc/W increases moves towards the Dirac point at E = 0, 
yielding an increase in Lc /W as a function of W for an 
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FIG. 4: The dependence of the conductivity on energy and 
length of a nanoribbon 272 atoms wide for V = 1.5. 



energy close to the Dirac point. In contrast, the contour 
plot for V = 2.2 shows that for all energies Lc/W de- 
creases with the width of the system. This implies that 
all states are exponentially localized for V — 2.2. The 
case of V = 2 is interesting since it corresponds to the 
crossover between the two behaviors. Indeed, for V = 2 
and close to the Dirac point, the contour lines of con- 
stant Lc/W are almost constant in energy, indicative of 
a critical behavior, in contrast to the V — 1.5 case, where 
the contour lines decrease in energy with increasing W, 
which is opposite to the V = 2.2 case. 

The most interesting consequence of the existence 
of the metallic-to-insulating transition for disorder 
strengths smaller than V ~ 2 is the possibility to have 
an on/off ratio as a function of density (gate voltage), 
which is infinite when taking disorder into account. In 
order to illustrate this point, we evaluated the average 
resistivity p = (R)y/3W/4L for V = 1.5 and plotted its 
inverse (c) in figure [5] as a function of energy and length 
of the nanoribbon for a fixed width of W=272 atoms. The 
numerical coefficient in determining p is the geometrical 
factor associated with the way we defined the atoms for 
our numerical implementation. Under our scheme, the 
total number of graphene atoms is given by W x L. We 
used (R) and not (G), because in most experiments p is 
measured as a function of density using a fixed current. 

For the maximum length (L=300), the conductivity 
vanishes at E = due to localization at the Dirac point. 
Moving away from the Dirac point the conductivity in- 
creases with energy due to the metallic behavior, be- 
fore reaching a maximum, which leads to a diverging 
fmax/%i n as a function of energy (density) for a given 
geometry and disorder strength. At even higher ener- 
gies, the conductivity decreases again due to localiza- 
tion. This behavior is consistent with some recent ex- 
periments, where a similar behavior has been observed 
by photo-emission, consistent with a metal to insulator 
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transition [3l[. In epitaxial graphene, a metal to insula- 
tor transition was also observed by molecular doping • 
Several authors have computed the inverse participation 
ratio to show the existence of a gap close to the Dirac 
point 33]. In a related work, the density of states was 
computed numerically and found to be consistent with 
a localization-delocalization transition close to the Dirac 
point (SSI- 

Summarizing, we have shown that there is a metallic- 
to-insulating transition in disordered graphene monolay- 
ers, which happens close to the Dirac point and leads 
to a window of energy, where a metallic behavior exists. 
This metallic behavior exists only for small enough dis- 
order (V < 2) and yields two mobility edges, one close 
to the Dirac point and another at a higher energy and 
dependent on the disorder strength. As a consequence, it 
opens the door for graphene based devices, which show 
an arbitrarily large on/off conductivity ratio. 

The author acknowledges financial support from 
NSERC, FQRNT, RQMP and INTRIQ. 
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